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Abstract 



The effective nonlinear response of films of random composites consisting 
of a binary composite with nonlinear particles randomly embedded in a lin- 
ear host is theoretically and numerically studied. A theoretical expression 
for the effective second harmonic generation susceptibility, incorporating the 
thickness of the film, is obtained by combining a modified effective-medium 
approximation with the general expression for the effective second harmonic 
generation susceptibility in a composite. The validity of the theoretical results 
is tested against results obtained by numerical simulations on random resistor 
networks. Numerical results are found to be well described by our theory. The 
result implies that the effective-medium approximation provides a convenient 
way for the estimation of the nonlinear response in films of random dielectrics. 
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I. INTRODUCTION 



Nonlinear response of random composite materials has been a subject of intensive re- 
search for many years [1]. The nonlinear response in such composites may be affected by 
several factors. For example, the response can be greatly enhanced by local field effects, and 
the percolation effect could change the character of the system [2-9]. In particular, enhance- 
ment of nonlinear optical response in random dielectrics may be useful for the design of novel 
switching devices in photonics and real-time coherent optical signal processors. A number 
of theoretical works have been developed specifically to weakly nonlinear composites [2-13]. 
An excellent review on the recent developments of the subject has recently been given by 
Shalaev [14]. 

While most of the work have focused on nonlinear effects such as the Kerr effect, Hui et al 
[15,16] have derived general expressions for the effective susceptibilities for second harmonic 
generation (SHG) and third harmonic generation (THG) in a binary composite of random 
dielectrics. In the dilute limit, the effective SHG and THG susceptibilities obtained in Refs. 
[15] and [16] reduce to earlier expressions found by Levy et al. [17]. For many applications 
and materials fabrication, samples in the form of thin films are essential. It is, therefore, 
useful to study the effective SHG in films of random dielectrics. In the present work, we 
consider films of random dielectrics consisting of spherical particles with nonlinear response 
embedded randomly in a linear host. We study the dependence of the effective SHG on the 
film thickness and the concentration of the nonlinear component. By invoking a modified 
effective-medium approximation (EMA) to incorporate the film thickness in the estimation 
of the local fields, together with our previous general expression for SHG susceptibility, the 
effective SHG susceptibility including the effect of film thickness can be calculated. The 
validity of our theory is tested against results obtained by numerical simulations on random 
nonlinear resistor networks of different thickness consisting of two different kinds of resistors. 
It is found that our EMA provides a good description of numerical results for the whole range 
of thicknesses corresponding to systems spanning from two dimensions to three dimensions. 

The plan of the paper is as follows. In Sec. II, we present our formalism on the effective 
SHG susceptibility of films of random dielectrics. Section III gives a description on the model 
of numerical simulations. Results of numerical simulations are compared with theoretical 
results in Sec. IV. 

II. EFFECTIVE MEDIUM APPROXIMATION FOR SHG IN FILMS 

We consider a macroscopically inhomogeneous medium consisting of a random mixture 
of two types of materials. The materials a and b have different macroscopic D-E relations. 
If we only include the quadratic nonlinearities, the general form of the D-E relationship at 
zero frequency would be 

D i = J2 e ii E i + dijkEjEk, i = x,y, z, (1) 

3 jk 

where Di (Ei) is the ith component of the displacement (electric) field D (E). For simplicity, 
we consider a composite in which a volume fraction p of nonlinear material a is mixed with a 
volume fraction 1—p of linear material b, i.e., (^(r) = d1 jk for regions occupied by material a 
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and dijk(r) = for regions occupied by material b. We also assume that the linear dielectric 
constants of both components are isotropic, i.e., = e a 5ij, with a — a (b) for regions 
occupied by material a (6). In general, when a monochromatic external field is applied, the 
nonlinearity of the component inside the composite will generate local potentials and fields 
at all harmonic frequencies. Here, we focus on the effective second harmonic generation 
(SHG) for a thin of thickness L. 

For a finite frequency external applied field of the form 

E (t) = E ^e-^ + c.c., (2) 

the effective SHG susceptibility can be extracted by considering the volume average of the 
displacement field at the frequency 2u in the inhomogeneous composite medium. A general 
expression for the effective SHG susceptibility has been derived in Ref. [15]. The result is 
that to first order in the SHG susceptibility dij k) the effective SHG susceptibility is given by 
[15] 

(d-2oj;ui,ij)ijk =< ^2uj-ildl mn K^ m jK^ nk >, (3) 

where < • • • > denotes a volume average over the inhomogeneous medium. The tensors K 
and K T are possible local field enhancement factors of the forms 

A 2W ;4X) = — , (4) 

and 



Eo,2ui;i 



= (5) 



with a similar expression for K^. nk (x). The tensor IfJ -(x), for example, is the induced mth 
Cartesian component of the electric field at position x and frequency u when a field Eq jU is 
applied in the j-th direction at the same frequency in a linear random medium of the same 
spatial inhomogeneous structure. Similar expression for the effective THG susceptibility can 
also be derived [16]. 

In a composite in which only the material a is nonlinear, the effective SHG susceptibility 
takes on the simple form 

^ rpa J« pa pa ^ 

, e _ ^ < ^2u; Q -2^,^A - ' / fi v 

U -2w;w,w — V p rp p K°) 

It is desirable to further develop a simple analytic approximation for the average of the 
local field factors in Eq.(6). One such approximation is the effective medium approximation 
(EMA). One possible effective medium approximation is to calculate the electric fields within 
each particle as if that particle is embedded in an effective medium with dielectric constant 
e e [16]. In this approximation, Eq.(6) becomes 

T)f e T)f e 
H e =nd a ( 2u> )( — ^ ) 2 (7) 
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where D is the spatial dimension and the subscripts uj and 2u specify the frequency at which 
the local field is evaluated. 

The effective linear dielectric constants and t\ u can be evaluated from the usual EMA 
for linear random composites by 

e a — e e 

a=a,b v / 

where p a is the volume fraction of the ath component and the summation is over the different 
components in the composite. For a thin film of thickness L, the behaviour of the system 
lies somewhere between D = 2 for small thickness and D = 3 for large thickness. To 
incorporate the effects of finite thicknesses, it is most convenient to convert the continuum 
description into a discrete version of the problem. In the discrete version, the system is 
considered to be a lattice consisting oi N x N x L sites, where L is the number of layers. For 
a D-dimensional hypercubic lattice, the coordination number z, i.e., the number of nearest 
neighbors, is related to D through D = z/2. For a single layer (L = 1), z — 4 and for a 
large number of layers [L — > oo), the system is three dimensional with z — 6. For a system 
of finite thickness L, the average coordination number becomes L-dependent and takes on 

z(L)= 2(3-1). (9) 

Incorporating Eq.(9) into Eqs.(7) and (8) by replacing the dimensionality D by z(L)/2 
constitute a set of equations for a modified EMA for the effective SHG susceptibility of films 
of random dielectrics. 



III. NUMERICAL SIMULATIONS USING RANDOM NONLINEAR RESISTOR 

NETWORK 

To test the validity of the modified EMA for SHG susceptibility in a film of random 
dielectrics, we follow the approach in Ref. [16] and perform numerical simulations using 
a random nonlinear resistor network. Such an approach has been applied previously with 
success to study linear optical properties and nonlinear properties such as the Kerr effect 
in random composites [18-21]. To model a system of finite thickness, two types of bonds 
are placed randomly onto a iV x N x L lattice, with L running from 1 to N modelling a 
2D system (L = 1) and 3D system (L = N), respectively. In the resistor network model, 
it is more convenient to define the model in terms of conductance of the bonds. With 
probability p, the bonds in the lattice are occupied by nonlinear components characterized 
by the two parameters: the linear conductance g a and the nonlinear susceptibility d, with 
the latter representing the response to two voltages at frequency uj. With probability (1 —p), 
the bonds are occupied by linear components characterized by a linear conductance g b . To 
model a mixture of metallic and insulating components at finite frequencies, we take g a and 
g b to be of the forms 

„. , 1 + iujRC -uj 2 LC . , 

0^ = R + M ' (10) 
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and 



g\u) = iujC. 



(11) 



The nonlinear susceptibility d of material a is assumed to be independent of frequency, for 
simplicity. This choice of g a and g b has been widely used in studying the frequency response 
in random composites [16,20,21], because the ratio g a /g b has the same form as the Drude 



dielectric function of a metal with the plasma frequence oo p = yl/LC and relaxation time 
r = L/R. In what follows, the parameters are chosen to be C — 1,L = 1, and R = 0.1, 
corresponding to a plasma frequency oj p = 1. 

The effective response of the random nonlinear resistor network can be studied by direct 
numerical simulations [16]. Using the Kirchhoff's law, the voltage at each node can be solved. 
Although the applied voltage imposed on the nonlinear network is taken to consist only of 
a component at frequency cu, the voltage at each node inside the network has components 
at frequencies uj and 2uj. So, solving the Kirchhoff's law leads to solving two equations 
J2 /H = and I^ 2 ^ — at each node simultaneously, where the summation is over all 
the bonds connected to the node under consideration and 1^ and /( 2w ) are the component 
of the current corresponding to the frequency u and 2u, respectively. These two equations 
at each node must be solved self-consistently for the voltages at frequencies uj and 2u>, and 
the effective SHG can be then extracted. In what follows, we use a lattice of size 20x20xL, 
with an applied voltage set to unity. By varying the layer thickness L, dimensional crossover 
of the effective SHG response from 2D to 3D can be studied. 



Numerical simulations allow for the extraction of the effective linear and nonlinear re- 
sponse of the inhomogeneous system for different concentrations of the nonlinear component 
and for different film thicknesses. Results from numerical simulations can then be compared 
with those calculated via the modified EMA so as to establish the validity of the approx- 
imation. Figure 1 and Figure 2 show the effective linear response g e ff and effective SHG 
response d e /d of the random nonlinear network as a function of frequency uj for three val- 
ues (p = 0.1,0.5,0.9) of the concentration p. In each panel of the figures, the symbols 
give the simulation results and the lines give the EMA results. The symbols represent an 
average over 50 independent runs corresponding to different configurations of the random 
network. The EMA results are obtained by using Eqs.(7)-(9), with the dielectric constants 
replaced by the conductances. Four different layer thicknesses are studied for each value of 
the concentration. 

In the dilute case (p = 0.1), the concentration is lower than the percolation threshold 
regardless of the thickness, i.e., p < p c for both p c (2D) and p c (3D) in the bond perco- 
lation problem. In this case, isolated nonlinear bonds (component a) are surrounded by 
the insulating bonds (component b). In the continuum case, the situation corresponds to 
one in which isolated disks (or spheres) of component a are embedded in a background of 
component b. The structure observed in g e ff corresponds to the so-called surface plasmon 
resonance structure. The surface plasmon resonance peak shifts to lower frequency as the 
layer thickness increases. The shift is related to the change in the resonance frequency from 




IV. RESULTS AND DISCUSSION 
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for an isolated sphere in 3D (large L) to uj p /\/2 for an isolated disk in 2D (L = 1). 
The real and imaginary parts of the nonlinear response normalized to the SHG coefficient 
of the nonlinear component d e j d show a corresponding resonance structure near the surface 
plasmon resonance frequency. The modified EM A captures the features reasonably well, 
with a better agreement for large values of L. It is reasonable in that the EMA, being a 
mean field approximation, is expected to work better in high spatial dimensions in which 
fluctuation effects are smaller. We also note that the simulation results indicate that the 
system picks up three dimensional character readily for not-so-large layer thicknesses. 

For p = 0.5, which is the percolation threshold for an infinite 2D system and higher 
than p c (3D), the numerical results for Re(g e ff) with L = 1 show a broad structure with 
the values at low frequencies tend to be enhanced. For L > 2, large Re(g e ff), relative 
to the dilute case, at low frequencies indicates that the composite behaves effectively as 
a conducting (metallic) system. The greater the thickness, the larger is Re(g e ff). It is 
because the value of p = 0.5 becomes increasingly larger than p c as L increases, and the 
system shows increasingly dominating metallic character. The effective SHG response shows 
a broader structure with an enhanced d e /d compared with the dilute concentration case. 
The enhancement in the SHG response comes about from the enhancement in the local 
field applying on the nonlinear component when the system has more better conducting 
nonlinear component. The modified EMA, again, gives results in reasonable agreement with 
simulation results, indicating that incorporating a thickness-dependent coordination number 
into EMA can be used to estimate the effective SHG response in random dielectric systems. 

For high concentrations (p = 0.9) of the nonlinear component, the systems shows metal- 
lic behavior in the linear response, as shown in Fig.l. The behavior of the effective SHG 
susceptibility basically follows that of the case of p = 0.5, with the appearance of a small res- 
onance structure at low frequency. The structure is due to the isolated insulating component 
in an otherwise metallic background. 

In both our numerical and EMA results, the SHG composite susceptibility is enhanced 
only relatively weakly over that of the pure SHG material. By contrast, such processes as 
Kerr nonlinearity are greatly enhanced in composite media [2-5,7-13,22,23]. We believe this 
weak enhancement may occure because Eq. (6) involves the linear fields at two different 
frequencies. Thus, the enhancement will be smaller than when all the field enhancement 
factors are resonant at the same frequency. The enhancement is further weakened by ran- 
domness when both components are present in high concentrations, which causes the field 
enhancement factors to be broadened over a range of frequencies. At very low concentra- 
tions, despite the two frequencies, other studies have shown that the SHG enhancement is 
still large. 

Also, our calculations assume that only the metallic component has a nonzero SHG 
susceptibility. We have done a few test EMA calculations in 2D, and find a somewhat 
larger enhancement, when only the nonmetallic component is nonlinear, with lm(d e ff/d) 
and Re(d e ff/d) reaching magnitudes as large as 4 or 5. 

Expression (6) for the effective SHG susceptibility treats the composite as a homogeneous 
medium with effective properties, as is reasonable if the inhomogeneities are small compared 
to a wavelength. Since the medium is effectively homogeneous, the SHG emission will be 
in the form of a collimated beam. Besides this beam, there is another SHG contribution, 
resulting from scattering from the inhomogeneities. This contribution, which has a broad 
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angular distribution, depends on particle size, unlike that calculated here. It has been seen 
in experiments on inhomogeneous metal-insulator films [24], and may be large near the 
percolation threshold. 

V. SUMMARY 

We have studied the dependence of the effective SHG response in a random dielectric 
system as a function of the concentration of the nonlinear constituent and the systemn 
thickness. A modified EMA, based on a previous expression of the effective SHG derived 
by the authors coupled with a thickness-dependent coordination number, is proposed for 
the effective SHG response. Numerical simulations on a random network consisting of a 
mixture of nonlinear and linear conductances allow extraction of the effective linear and SHG 
response. It is found that the modified EMA gives results that are in reasonable agreement 
with the numerical simulation results. The system picks up three dimensional character 
for moderate film thicknesses. The results show that the modified EMA can be used as an 
effective tool for estimating the effective SHG response in random dielectric systems. The 
approximation can be applied to other systems consisting of nonlinear components. For 
example, the nonlinear SHG response may be associated with the poorer conductor of the 
two constituents. In this case, the contrast in the linear dielectric constants between the 
two constituents may produce an enhanced local field in the poorer conducting component 
and thus drive a greater SHG enhancement. Indeed, some preliminary calculations, using 
the EMA, do give rise to this increased enhancement. 
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FIGURES 



FIG. 1. The effective linear response g e ff as a function of uo. The three panels correspond to p 
= 0.1, 0.5, and 0.9. In each panel, results for thickness L = 1, 2, 4, and 6 are shown. The symbols 
are results of numerical simulations and the curves are theoretical results based on the EMA. 

FIG. 2. The effective SHG response as a function of u>. The three panels correspond to p = 
0.1, 0.5, and 0.9. In each panel, results for thickness L = 1, 2, 4, and 6 are shown. The symbols 
are results of numerical simulations and the curves are theoretical results based on the EMA. 
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